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Abstract 

Using an internal-variable formalism as a starting point, we describe the vis- 
coelastic extension of a previously-developed viscoplasticity formulation of the com- 
plete potential structure type. It is mainly motivated by experimental evidence for 
the presence of rate/time effects in the so-called quasilinear, reversible, material re- 
sponse range. Several possible generalizations are described, in the general format 
of hereditary-integral representations for non-equilibrium, stress-type, state vari- 
ables, both for isotropic as well as anisotropic materials. In particular, thorough 
discussions are given on the important issues of thermodynamic admissibility re- 
quirements for such general descriptions, resulting in a set of explicit mathematical 
constraints on the associated kernel (relaxation and creep compliance) functions. In 
addition, a number of explicit, integrated forms are derived, under stress and strain 
control to facilitate the parametric and qualitative response characteristic studies 
reported here, as well as to help identify critical factors in the actual experimental 
characterizations from test data that will be reported in Part II. 

Keywords: viscoelasticity, hereditary behavior, TIMETAL 21S, nonisothermal, 
deformation, multiaxial, thermodynamics 
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1 Introduction 


A number of advanced material systems (for example metallic, polymer and ceramic 
based systems) are currently being researched and evaluated for high temperature air 
frame and propulsion system applications. As a result, numerous computational method- 
ologies for predicting both deformation and life for these classes of materials are under 
development. An integral part of these methodologies is an accurate and computation- 
ally efficient constitutive model for the matrix constituent in such systems. Furthermore, 
because of the proposed elevated operation temperatures for which these systems are 
designed, the required constitutive models must account for both time-dependent and 
time-independent deformations. For example considering that most aerospace engine 
designs are typically limited to the quasilinear stress and strain regimes, the reversible 
time-dependent response component becomes dominate in comparison to the irreversible 
component. Alternatively, one can envision another extreme case (e.g., in polymer and 
rubber based systems under varying temperatures) in which a purely reversible viscous 
response is present. And lastly, an obvious natural extension for general applicability 
is the middle ground in which a combined reversible and irreversible representation is 
required. 

To accomplish this we will extend a previously developed, complete potential based, 
framework [1],[2] utilizing internal state variables [[3], [4]] which was put forth for the 
derivation of time-independent reversible and time-dependent irreversible constitutive 
equations. This framework, and consequently the resulting constitutive model, is termed 
complete because the existence of the total (integrated) form of the Gibb’s complemen- 
tary free energy and complementary dissipation potentials are assumed a priori. In 
outline form, expressions for the Gibb’s thermodynamic and the complementary dissipa- 
tion potential functions are assumed in terms of a number of state and internal variables 
characterizing the changing internal structure of the material. For instance given the 
Gibb’s potential in the following form 


i = ( 1 ) 

and assuming a priori that the inelastic strain is an independent parameter (and 
not an internal state variable), for example 


<?kk 


$ = E(<r ijt T) - < 7 * 4 + H K,T) - -fu(T - T 0 ), 

an expression for the total strain rate can be obtained by differentiating, that is, 


_ d ( ~d$ _ . 

<J ' " dS da fj * “ + ij + 


ddijdT 


Sij ( du> . _Y 

3 \ dT )\ 


as well as the rate of change of the conjugate internal variables (A^), 

d —(94> 


( 2 ) 

( 3 ) 

( 4 ) 
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where 


( 5 ) 


_ -a»* _ 

ijrs daijda,.. 


and 


Qil ~~ 


d<Jxj da r3 


da^dou, da^doUf 

are the external and internal compliance operators, respectively, and 


B, = 


_ -fr H (q 7 ,T) 


( 6 ) 


(7) 


da^dT da^dT 

is the change in the conjugate internal variable (.A^) with temperature. Note the three 
terms in eqn. (3) may then be identified from left to right as the elastic (time-independent 
reversible) , inelastic (irreversible) , and thermal expansion components of the total strain 
rate, respectively. Thus, 


e-- = g -)- e'f. -|- cF. 

c ij ' tj • t] 


where 


and 


with 


and 


— Cijrs&rs 


= MijT 


-&<*> 
ij ~ dcTijdr 


Mu = 


-frEfa) Sij 

d^&f + 3 


MX = W + — (T — T 0 ) 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


denoting the instantaneous coefficient of thermal expansion, and e£ (the inelastic strain 
rate) is defined separately using the concept of a complementary dissipation potential 
Q(a i j,a 7 ,T). Whereby, given an 


Q. = fi(<ty > a r ,T) 

and using the Clausius-Duhem inequality[5]; the flow law becomes 

an 


cf = — 

a* 


(13) 


(14) 


V 


and the evolutionary laws for the thermodynamic conjugate internal state variables are: 
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an 

*^ 7 da-. 

Utilizing eqn. (4) the internal constitutive rate equations for the internal state variables 
are obtained, 

= Lyi [A< -W] ( 16 ) 

where 



Thus, eqns. (14) and (15) represent the flow and evolutionary laws, for an assumed 
0 = n(<Tij,<Xy,T), and eqn. (16) the internal constitutive rate equations, given a Gibb’s 
potential $, wherein both potentials are directly linked through the internal state vari- 
ables Cfcy. The specific forms selected thus far for both the Gibb’s and complementary 
dissipation potentials resulted in a fully associative, multiaxial, nonisothermal, unified 
GVIPS (Generalized Viscoplasticity with Potential Structure) model with nonlinear 
kinematic hardening [6], [7]. 

During the detailed experimental program to specify the required material functions 
and characterize the associated material parameters for TIMETAL 21S 1 , an advanced 
titanium-based matrix commonly used in titanium matrix composites (TMCs), it was 
discovered that TIMETAL 21S exhibited both a time and temperature dependent re- 
versible (linear viscoelastic) and irreversible (viscoplastic) domain. These reversible and 
irreversible domains are posited to be delineated by temperature dependent threshold 
surfaces, as illustrated in Figure 1 for low, mid and high temperature regimes. In partic- 
ular, as opposed to the purely time independent (elastic) behavior within the inner most 
threshold surface of Fig. lb, two additional time dependent (reversible and irreversible 
mechanisms) will be activated upon traversing the two threshold surfaces shown. As our 
previous work has focused on the specification and characterization of the irreversible 
domain (i.e., eqns. (13)-(17)) the primary objective of the present study will be to con- 
struct, within the context of a complete potential structure, a multiaxial, nonisothermal 
viscoelastic model to describe the reversible strain component, see eqn. (9), of the total 
strain decomposition. 

Clearly, this extension represents a very difficult problem [[8], [9]] in view of the multi- 
tude of choices available through the general functional forms $ and D, and the complex 
interaction between the newly desired reversible viscous response component and the 
previous elasto-viscoplastic contributions. In particular, this involves the selection of 
the type of additional internal state variables and/or parameters accounting for the new 
mechanism, and their corresponding implied partitioning format of the stress and strain 
variables (further elaborated on in the following perspective section), within the context 

1 TIMETAL 21S is a registered trademark of TIMET, Titanium Metals Corporation, Toronto, OH. 
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of miiltiaxial stress states and in conformity with any thermodynamic admissability 
requirements (see Section 5 and references [[8], [9]] for many pertinent controversial 
issues). 

The specific model construction has been motivated by experimental observations as 
discussed in detail in Part II [10] of this report. As an example see Fig. 2 which illus- 
trates the rate dependence of the uniaxially obtained Young’s Modulus as a function of 
temperature. Clearly at elevated temperature the Young’s Modulus becomes increasingly 
rate dependent, thus indicating a need to include some type of viscoelastic influence into 
the deformation model at these temperatures. Another important observation, as shown 
in Fig. 2 of Part II [10], is the purely transient time dependent behavior of TIMETAL 
21S with no steady state behavior being exhibited within the reversible strain domain 
of this material. 

Other observations of a similar nature as above, i.e., marked rate and time-dependency 
in the so-called quasilinear range (at relatively low stress levels) have also been made by 
other investigators, motivating a number of extensions in constitutive theories. For ex- 
ample, these include distinctive static and dynamic moduli, reflecting the load-frequency- 
dependent behavior [11]-[12], and the use of an ’’initial” primary creep response compo- 
nent as extension of the well-known model of viscoplasticity theory based on overstress[13] . 

An outline of the remainder of the paper is as follows: In section 2 we discuss the 
significance of deciding a priori how one will partition the required stress and strain 
components and how this selection will influence our overall objective. In section 3 a 
specific multiaxial, nonisothermal, linear viscoelastic model of the internal-state- variable 
type is developed and then reduced and analyzed in section 4. In sections 5 and 6 
numerous generalizations pertaining to both the theoretical foundation and numerical 
implementation are discussed. Lastly in section 7 a parametric study is conducted using 
the simplification (section 4) of the specific model developed in section 3 to illustrate the 
importance of a key assumption (i.e., equality of the Poisson’s ratio’s) and to identify 
important factors in the actual characterization of the specific-model put forth. 


2 Perspective 

Traditionally, for the small-deformation problems considered here, the total strain has 
been partitioned into an elastic (reversible), an inelastic (irreversible) and thermal (re- 
versible) strain component, that is: 


c *j c t j ' ij ' 


(18) 


where the inelastic strain corresponds to such physical phenomena as tune-independent 
plastic strain or time-dependent viscoplastic strain (sometimes referred to as creep strain). 
Whether a material’s phenomenological behavior is simple or complex, a mechanical 
model is often of great help in the visualization of that behavior. Here we want to 
discuss four such one-dimensional mechanical representations so as to put into perspective 
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previous work and motivate the construction of the present viscoelastic model. To this 
end, as shown in Fig. 3, a viscous element is shown pictorially by a dashpot whereas a 
linear-elastic element is pictured as a spring. For each element of the mechanical model, 
change in length represents strain in the material and force represents stress. 

Let us begin with a uniaxial mechanical model that is representative of numerous in- 
ternal state variable viscoplastic models available in the literature (e.g.,[4]), in particular 
the recent unified GVIPS model put forth by the authors[[6], [7]]. Fig. 3a presents a 
generalized nonlinear four element model (a spring in series with a nonlinear Maxwell 
element in parallel with a nonlinear dashpot) in which the strain is partitioned into an 
elastic and inelastic strain component; wherein implicit in the inelastic strain component 
is a partitioning of the stress into an internal stress (a) and effective stress (a — a ) com- 
ponent. The resulting system of differential equations (total strain, flow, and evolution 
laws) are also shown in Fig. 3a, where it is clear that the flow and evolutionary laws are 
identical in structure (provided that the material parameters (e.g., Ei, T72) are appropri- 
ately modified to become state dependent) to the GVIPS form, e.g., in the evolutionary 
law a competitive mechanism between hardening and thermal recovery is present. 

Conversely, if we examine the classical four element (non)linear viscoelastic model 
[[14], [11]], as illustrated in Fig. 3b, which is capable of producing a similar elastic, primary 
and secondary creep history, we see that the resulting form of the flow and evolution 
laws are quite different. This difference is predicated upon the further partitioning of 
the inelastic strain component into a primary and secondary component, achieved by 
placing a Kelvin model in series with a dashpot, thus resulting in a more complex flow 
law and restrictive evolution law. Comparing the resulting equations of Fig. 3b with 
that of Fig. 3a, it is immediately apparent that one should not assume (as is often 
done) that commonly employed unified viscoplastic formulations (of which Fig. 3a is 
representative) are merely extensions of this classical four element configuration to the 
irreversible domain. 

As the emphasis in this study will be on developing a transient viscoelastic model, 
let us compare two three element configurations that give the same "classical” 2 uniaxial 
behavior in both creep and relaxation. The first can be immediately obtained from a 
truncation of the four element model presented in Fig. 3b. This truncation is obtained 
by setting T74 to infinity and results, as before, in a partitioning of the applied strain into 
an elastic and primary creep regime. Alternatively, one can construct a complementary 
configuration in which the stress (instead of the strain) is partitioned into an elastic and 
viscous part. This model is known as the standard linear (n = 1) solid model[ll] and 
is shown in Fig. 4a along with the pertinent one-dimensional governing equations. The 
one dimensional solution to the governing systems of equations for Figs. 3b and 4a, 
for the special case of linear (n = 1) viscoelasticity, is illustrated in Fig. 4b for both 
models and will be subsequently used to validate the one dimensional simplification of 

2 By classical we mean, in conformity with the classical rheological model we interpret the uniaxial 
results to imply complete neglection of the strain component interactions (i.e., no Poisson’s effect and a 
truly one-dimensional response) 
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the proposed multiaxial theory. Clearly, both approaches (i.e., partitioning of the stress 
or strain) are capable of producing the desired uniaxial response (see the table in Fig. 
4b); however, with an eye toward integrating the resulting viscoelastic theory with the 
author’s previously proposed GVIPS model, the partitioning of the stress (standard 
linear model) approach is preferred. 

This preference is obvious when one recalls the four element model of Fig. 3a, the 
decomposition of the total strain into reversible and irreversible components, and our 
desire to maintain the current GVIPS representation for the irreversible component of 
strain. Consequently, we will merely replace the single spring element in Fig. 3a with the 
standard linear model of Fig. 4a, as illustrated in Fig. 5. Mathematically this is accom- 
plished by employing an additive decomposition for the two underlying thermodynamic 
potentials (i.e., complementary energy/ Gibb’s (see eqn. 1) and dissipation functions (see 
eqn.. 13) into reversible and irreversible parts, that is $ = $r + $ir and fi = Qr 4- fiiji; 
where $ IR = H K,T) and Q IR = f l(ofc, a^,T) are defined as before in eqns. 

(2) and (13). Henceforth, we will concentrate only on the construction of the multiax- 
ial, nonisothermal, linear viscoelastic representation of the now (tune-dependent) viscous 
reversible strain component, i.e., and Hr . 

3 Multiaxial Nonisothermal Theory 

Here we will restrict our discussion to the reversible viscous part, where the corresponding 
functions <3 >r and Q.r are assumed a priori to be in conformity with Fig. 4a. Further- 
more, for conciseness the discussion is limited to a case involving small deformations (in 
which the initial state is assumed to be stress free) and a Unear viscous element (which 
is in accordance with the experimental evidence described in Part II). A Cartesian coor- 
dinate reference frame and index notation are utilized (repeated Roman subscripts imply 
summation). 

The reversible contribution to the Gibb’s potential takes the following partitioned 
form 


Otj, T) = ~ 


m — — 
T8 


' kk 


+ o] 


kk 


w(T-T 0 ) (19) 


where the first term represents the stored energy in the spring element, the second, 
the stored energy in the Maxwell element, and the third, the energy dissipated by the 
dashpot within the Maxwell element. The equations of state for the conjugate strains 
(i.e., strains in the spring and Maxwell elements) can be obtained by differentiating 
with respect to the associated variables, and < 7 ™ , respectively; i.e., e^- ~daf^ an d 

e ij = • 111 rate form > 


e*. = _ ( — g) = E~ 


datj 


-1 Z.9 
ijT6°ra 


+ 0%T 


( 20 ) 
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( 21 ) 


;m _ d f 1 m 4 a Mrfi 

ij Jt K ~d^r ] ~ Mi ’ r * rg + e * + 0i > 1 




pM _ 

0 ij ~ 


* , ^y 

&f " + ¥"“ 


-3M U !, 

ST 


^ra “I" o ^tan 


where denotes the instantaneous coefficient of thermal expansion (see eqn. (12)), and 
e?. (the dashpot strain rate) is defined separately using the concept of a complementary 
viscoelastic dissipation potential f2«(cr™, T). 

Given 

n* = (m) 

the flow law for the dashpot becomes 


■d 9 CIr ! m 

a n m 7 lijrs°rs 
u u ij 


Thus substituting eqn. (25) into eqn. (21) and making use of the fact that 

= o„ - < 7 ?. (26) 

due to the partitioning of the stress, eqn. (21) becomes, 

e£ = M£,(& rs - a s rs ) + r? t ^ a (cr„ - O + 6%T (27) 

As a kinema tic compatibility constraint, the total strain is equivalent to the strain in the 
spring and in the Maxwell element, i.e., 

ey = = e«, (28) 

the evolution of the internal stress(a* J ) can be derived. 

K, = - O + e v t} (29) 

where 


EM uij = (E„;, + M t -'} 


+ ^(<7.. - O (31) 

Finally, using eqns. (20) and (29) we obtain the desired multiaxial nonisothermal vis- 
coelastic stress-strain expression, 
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(32) 


&ij = Eijh [EM klmn { M m n ra &ra + Vvmrs^ra - &**) + W}] + 

•with the evolution of internal stress given by eqn. (29). 

Now assuming an isotropic material, the general form for the isotropic material 
tensor (in terms of the bulk modulus K and shear modulus G) is as follows: 


C ijkl = (. K c - | G C )M« + 2 Gci^ji + S it 6 jk ) (33) 

where the subscript on the bulk and shear modulus (i.e., C ) will always be indicative of 
the particular material tensor being replaced (e.g., Cym )- Similarly, one can write the 
general isotropic stress-strain relations as follows: 


Oij = Ke kk Sij + 2 G(eij — -e kk 6ij) (34) 

e H = ^ kk 6ij + ~ !»**$&■) ( 35 ) 


Utilizing these above relations and performing a number of algebraic manipulations on 
eqn. (32) the desired isotropic, multiaxial, nonisothermal, viscoelastic stress-strain ex- 
pressions can be obtained, that is 


J^E^kk^ij + 2Ge {&ij — ~£kkfiij) 


Kem 

Km Gm 


Gem 1 i. c . Gem ^ , 

CTkkOij + — — + 


Gm 


(36) 


1 Kem Gem], ju , Gem , n . 

X -T? 77" ~ °rr ) b H + ~n Wl + 


K, 


Gr, \ 


S - %] ^ - O'* <««+ 

- a Sf<, + \ ft - 1] ^ + 


where the notation for the derivative with respect to temperature of the bulk and shear 
modulus is defined as 


Ke, m — 


K%, m CO 
SKb, m(T) 

err 


(37) 


r T - 

^E, M — 


^E, M CO 
m(T) 
dT 


(38) 


such that if K and G are temperature independent, i.e., ^ = 0, K^,Kjj and G£, Gjj 
are infinity and thus the temperature rate term ( T ) will vanish. 
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4 Subspace Simplification 

If we consider a triaxial (principle) state of stress and isothermal conditions (T=0), eqns. 
(36) can be transformed into the following system 


Y 

X + AX 

with ( ) denoting the time derivative and 



DX 


AL + AL 


"Jl \ ; 

f a U 

a 22 I > L = | 

&22 

<^33 / 

\ <^33 


D = 


3 Kg + Ge 
9G b K e 


1 

(3X g -2G g ) 
2(3 K e + G e ) 
(3iC g -2G g ) 
2(3 K b + G e ) 


(3K g -2G g ) ( 3 K e - 2G e ) - 

2(3 K b + G e ) 2(3Jf g +G g ) 

, (3Jf g -2G g ) 

1 2(3JC g +G g ) 

( 3 K b - 2G e ) -i 

2{3Ke+Gjb) 



’ b 

a 

a 


’ b 

a 

a 

A = 

a 

b 

a 

, A = 

a 

b 

a 


a 

a 

b _ 


a 

a 

b m 


where 


(39) 


(40) 


(41) 


(42) 


1" Kem Gem 

, b=- 

K E m ^ Gem ] 

CD 


CD 

1 

!*? 

Kem Gem 

, b=- 

Kem 0 Gem 

k m g m 

’ 3 

. Km Gm . 


(43) 

(44) 


Now as the matrix A is real, symmetric and positive definite, then we know that 
1) its eigenvalues (Ai, A 2 , ...) are all real and positive, 2) its eigenvectors can be chosen 
orthonormal and 3) A = NAN -1 where A is a diagonal matrix containing the eigenvalues 
of A. Consequently, if we let 


Z = N -1 X, Z = N -1 X (45) 

eqn. (39) can be rewritten as 

Z + AZ = N -1 (AL + AL) (46) 


and solving for Z we obtain 


Z = Ce -At + Z por * c “ ior 


(47) 
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where the form of Z port ‘ cuiar will depend upon the applied loading vector L. Finally, 
transforming back to the real space X, we obtain the desired general solution 

X = NCe _At + X parti ™ lar (48) 


with 

Y = DNCe At + DX port<cu,or 

4.1 Creep History 

As an example, let us examine the case when the applied load is linear in time, i.e., 

L = 0t + L* 

1 = 0 

wherein the (3x1) vectors 0 and L* are prescribed time invariant constants; thus eqn. 
(46) maybe rewritten as 

Z + AZ = N -1 [A0t + A0 + AL*]. (49) 

Solving for the particular solution of eqn. (49) we obtain 

TP = A" 1 [(tl - A" 1 ) N~'A0 + N-\A0 + AL*)] (50) 

and transforming back to X we obtain 

X p = 0t- A~ l \0 - A0\ + L* (51) 

so that in this case _ 

X = Ne- At C + 0t- A _1 [/3 - A0] + L* (52) 

where the coefficient vector C must still be solved for, given a set of initial conditions. 

For instance assuming that we start at zero stress and apply a constant load rate and 
then at time equal t 0 hold the load constant (i.e., a typical stress controlled creep test) 
we would obtain the following, that is for: L* = 0 and X = 0, 0 ^ 0 at t = 0 we get 



C = N -1 A _1 [0 - A0] 

(53) 

and 

X = (Ne _At N _1 - 1) A" 1 [0 -A0]+0t- 0 < t < t 0 

(54) 

and for t>t o ,(0 = O and L* ^ 0) we get 



C = [Ne" At0 ] 1 [X^ - L*] 

(55) 

and 

X = Ne” At { [Ne- At °] _1 [X* - L*] } + L*; t h to 

(56) 


11 



where 

with 


X*o = (Ne~ At °N -1 - 1) A" 1 [/3 - A/3] + /3to 
Y = DX 


(57) 


as before. Or alternatively, we could assume the classic creep case involving an instanta- 
neous load-up to a constant stress value L* ^ 0 at to = 0 and therefore after taking the 
appropriate limits as to approaches zero (i.e., lim X^ = A/3 ) eqn. (55) becomes 


find eqn. (56) 


and if /3*= L* : 


C = N _1 [A/3*-L*] 


X = Ne -At N -1 [A/3*— L*] + L* 


X = {Ne -A *N -1 [A - 1] 4- I}L* 


(58) 


(59) 


(60) 


4.1.1 Uniaxial Simplification 

Now if we further restrict ourselves to a purely uniaxial creep loading (e.g., L\ = <7n = a* 
& &n = Pi = <722 = <733 = CT22 = <733 = 0) case and make use of the symmetry resulting 

from the previous assumption of an isotropic material (i.e., = <733, 0^22 = ^33)) the 

various matrices and vectors take on the following reduced forms: 


A = 
A -1 = 


b 2a 

1 

a a + 6 

“3 


K 

K 


• + 2 G* 2 K* - 2G* 1 
:* - G* 2K* + G* \ 


A = 


K* 0 
0 G* 


3G*K* 

N = 


G* + 2K* -2 (K* — G*) 

- (K* - G*) 2G* + K* 


1 2 

? 3 

75 7s 


N _1 = 


-v/lO 


-1 _ 2 

75 yS 

"75 75 


b 2d 1 _ 1 [ K* + 2G* 2K* - 2G* ] 
d d + b \ ~ 3 [ K* - G* 2 K* + G* \ 


where 


K* = 


K 


EM 


K n 


G* = 


Gem 
G„ ’ 


K* 


l tj 

K, 


EM 


G* = 


G 


M 


EM 

y 

T M 


(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 
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So that 


Ne -A 4 N _1 = \ 

o 


e~ KH + 2e~ G,t 2 (e'** 4 - e" 6 * 4 ) ' 
e -K‘t_ e -G‘t 2e~ K ' t + e -G * 4 


( 67 ) 


a* 


A 1 [/3 A/3] 9Q* K * to 

and in the classic creep case (i.e., eqn. (60)) 


[1 - (K* + 2 G*)] (2K* + G*) + 2(A* - G*)(tf* - G*) 
- [1 - (> + 2G*)] (A* - G*) - (A* + 2G*)(K* - G*) 


x =f 


3 + (K* - l)e~ K,t - 2(1 - G*)e _G 4 ) 
{K* - l)e _ir * 4 + (1 - G*)e -G * 4 
(K*- l)e-^ 4 + (l-G)e- G * 4 


and 


Es 



( 68 ) 


(69) 


(70) 


where 


K* = E*v, G' = E*v\ K* = E*u, G* = E*v* 


(71) 


„ _ EM _ / 

' EmEs \ 1 

EM _ 

E s 

V ' 

k Em + Es) f] 

Em 

Em + Es 

(1 ~ 2» '„) 

(1 — 2uem) ’ 

(1 + Vem ) ’ 

E m v s + E s vm 
vem — t, , ei 

&S + 


(1 - 2v m ) _ 

(1 + v M ) 


P = 

(1 — 2vem j ’ 

(1 + vem) 



(72) 

(73) 

(74) 


and D was taken to be 



(75) 


Now in order to recover the ” classical” one-dimensional linear viscoelastic solution 
for an instantaneously (step-function) applied creep load, see Fig.4b, all Poisson’s ratios 
must be taken to be zero (i.e., v s = v m = v n = 0, leading to a=0 and K* = G* = 
ft* — G* — E*). Consequently, eqn. (70) reduces directly to the simple expression 

given in Fig. 4b, that is 



1 - 


(E* - l)e _E * 4 
0 
0 


(76) 
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where E * (i.e., £) and E* are defined in eqn.(72). Alternatively, we can rewrite eqn. 
(70) in terms of individual Poisson’s ratios and element stiffnesses, to obtain the initial 
or instantaneous response (that is at t=0): 


Y = C o 


(1 + a) - (1 + 2u m )(u m + otv,) - 2i/ a [(v m + ai/,) - i/ m (l + a)] 
(u m + au s )(l + 2u m u a ) + (1 + a)(i /,( - 1) - i/ m ] 

(u m + ai/,)(l + 2v m v s ) + (1 + oc)[ug(i/ m - 1) - v m ] 


where 

c Zl 

0 E s [(1 + a) - 2(u m + au,)\ [(1 + a) + ( v m + ai/,)] 
Then if we assume v t ^ u m , but a = 1 then, 


(77) 


Y = 


(Y 

2E S [(2 - (u m + v a ) - (i v m + v a y] 


2-{v m + v a )-(2ij + 2vl) 

-i/,(l + i/ m )(l - 2u m ) + -i/ m (l + i/,)(l - 2u a ) 
-u a (l + I/ m )(l - 2u m ) + -I'mil + ^a)( 1 - 2l/,) 

(78) 


or if we take v a = v m = u, but a / 1 then, 


E a ( 1+a) ] 


1 

—v 

—v 


a* 

E a + E m 


1 ' 
-v > 
—u 


(79) 


Consequently, the above expressions provide us with an alternative means (instead of 
direct experimental measurement) to determine the proper assumption for the Poisson’s 
ratios, given the experimentally measurable dynamic modulus. It is interesting to note 
that if the assumption of constant and equal Poisson’s ratios is accurate then eqn. (79) 
proves that the dynamic modulus is equal to the summation of E a and E m , i.e., Eq = 
E a + E m . 


4.2 Relaxation History 

Alternatively, considering a triaxial (principle) state of strain and isothermal conditions, 
a similar set of expressions as those shown in eqns. (39)-(60) can be obtained. That is, 


where 


X = D -1 Y 

L + BL = BY + BD _1 Y 


B = A -1 A 
B = A -1 D -1 
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and A b = NbBN b with N B being the eigenvector matrix of B . Then assuming the 
applied strain field is linear in time, i.e., 


Y = C* + Y* 

Y = C 

we obtain the counterparts to eqns. (54), (56), and (60), respectively: 

L = (N B e _Ast Ng 1 - I) B- 1 [D" X C - BC] + D" 1 #; 0 < t < to (81) 

L = N B e- Aflt { [N B e- A * t0 ] _1 [L* - D -1 Y*] j +D -1 Y*; t h to (82) 

and 

L = {N B e' Ast N^ [B - D' 1 ] + D 1 } Y*, if to = 0 (83) 

The above general triaxial forms will prove very useful in conducting parametric 
studies under multiaxial states of stress and in providing insight into identifying the 
most discriminating experimental validation test conditions. In fact, all the results to be 
given later are generated using these principle field expressions. 


5 General Characteristics and Equivalent Hereditary 
Integral Representation 

Largely motivated by the experimental observations alluded to earlier (see Fig. 2 and Part 
II) we have elected to focus our attention on a single internal dissipative mode in the above 
derivations. However, the underlying thermodynamic formalism [[5], [9]] is sufficiently 
general to provide for a number of possible extensions. Several such generalizations are 
described below, together with the exploitation of their ” equivalent” integral forms within 
the modem theories of materials with memory; i.e., hereditary integral representations, 
e.g. [9]. In particular , a difficult question concerns here the proper conditions to be 
imposed on the kernels associated with general hereditary integrals, and in this connection 
the distinct advantages offered by the present formulation will become apparent. 

5.1 Limiting and Equilibrium States 

As a simple rheological interpretation the assumed structure of the (complementary) 
free energy, eqn.(19), leads to a local additive decomposition of the stress tensor into 
an equilibrium stress, u s , and a non-equilibrium stress, cr m , see chapter 2 of [15] 
and Kremple formulation in [13]. The non-equilibrium stress is governed by a linear, 
dissipative, evolution equation exhibiting the well known fading memory phenomenon.. 
This is characterized by a viscous-like environment in which microstructural changes 
take place within a temperature dependent quasi-linear regime, as illustrated in Fig. 1. 
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In particular, for arbitrary deformations, from eqns. (26) and (29), the condition for 
thermodynamic equilibrium is when 


<r m |t-»oo= 0 (84) 

This in turn results in ifj |t-»oo= 0 and hence the vanishing of the dissipation at equilib- 
rium, i.e., Qr | t ->oo= 0- This is in conformity with the thermoelastic (time independent 
reversible) limit put forth by Truesdell and Noll [5] wherein the thermodynamic process 
becomes reversible, and the material reacts fully thermoelastically. This limit is depicted 
schematically in Fig. 6 by the solid line, and the following one-dimensional expression , 
that is , 


= £o‘<rr;(T) + ££' <<r - <r„ (T)> (85) 

where Eq is the initial (instantaneous or dynamic) stiffness modulus, is the stiffness 
modulus at infinite time (E s ), Oti is the upper stress limit delineating time independent 
reversible behavior from time dependent, and () denote Macauley bracket. 


5.2 Model Generalizations 

For simplicity, and without loss of generality, only the isothermal conditions are consid- 
ered in the following discussions. Furthermore, for the case of isotropy (see also Remark 
2) we assume all the fourth-order tensors of moduli E, M, t) to be coaxial. This, to- 
gether with the linear structure of eqns. (27) and (29), enables us to express the internal 
variable, cr m , in terms of a simple convolution integral: 

e “ ^j-[ rf*'E~ l or]ds 

as 

or 

<r m = I e ' — fa*e] ds (87) 

Jo 

where EM = 7*7 = * 7 * , M = ^17 = rj*are assumed with r and p as the positive relaxation 
times. 

With this convolution representation many generalizations are possible once an 
appropriate choice of its kernel is made. For example, several relaxation times may be 
introduced by replacing the single exponential term in eqn. (86) with the discrete- 
relaxation-spectrum kernel, 

N 

K(t) = 'Z, er ‘ ( 88 > 

j = 1 

The most general form, that is a continuous spectrum of relaxation times could al- 
ways be used [[12], [8]] . Similarly <r m , in eqn. (86), maybe defined by convolution 
forms involving fractional derivatives [[16], [17], [9]], or any power type of kernels, 



( 86 ) 
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as in [[12], [9], [18]]. Furthermore, it has been demonstrated by Roscoe[19] that, as the 
number of terms increases, the discrete-spectrum model can approximate any given re- 
laxation/creep function with any desired accuracy, thus obviating the need to further 
study more complex rheologic models. In addition, very efficient general methods are 
presently available [e.g., [12] and [20]] to determine simplified continuous spectra 
(i.e., with small order or few terms that can very well approximate fairly complicated 
formulae fitted to the experimental data (e.g., creep law’s of the power or log-power 
type) Hence, from both the physical (mechanical-analogs) as well as the computational 
(i.e., recursive formulas with minimal storage requirements for numerical integration; 
see section 6 for details) standpoints the differential models (and their correspond- 
ing exponential-integral- forms, as in eqns. (86-88)), are certainly more appealing for 
applications.. 

Remark 1 As a consequence of the assumed coaxiality of material-moduli tensors E, 
M, and tj, the equivalent viscoelastic Poisson's ratio is rendered constant (i.e., time- 
independent), thus leading to a neat separation of the characteristic relaxation/retardation 
times p and t (see eqns. 86 and 87). Despite the lack of extensive experiments on 
this specific aspect , there is indication supporting a very small (if any) time-dependent 
change in Poisson's ratio . In addition , several other arguments in its favor have been 
also made, e.g., in the solution of boundary-value problems [8]. On the other hand, 
directly assuming a time-invariant Poisson's ratio is fully consistent with the important 
theoretical requirement of "damping-extent ordering"; i.e., damping in pure shear is 
larger than in extension, which in turn is larger than damping in bulk, a conclusion that 
can be reached without any consideration of Poisson's ratio. This is a direct consequence 
of the "compliance" nature of Poisson's ratio as a material parameter ; e.g., implying 
that the lateral contraction must lag behind the imposed sinusoidally oscillating strain in 
experiments involving harmonic excitations. This is equivalent to the conclusion that the 
imaginary part of the complex, frequency-dependent , Poisson's ratio must be negative; 
an issue that has resulted in much controversy in the earlier literature on this subject, 
e.g., see[[12], [21]] 

Remark 2 As a corollary to the above remark, one can also introduce material anisotropy; 
i.e., by simply utilizing the same anisotropic forms for the three tensors E, M, and rj. 
All other equations, particularly eqns. (86) and (87) remain unchanged. For example, 
with transverse isotropy, we have (see [22]) , where 

E = 2G,P -aR + 2(Gj-G t )Q + (A + | + 1 G t )6 0 6 + [3a + 4 (G, - G t ) + P]D 0 D (89) 
where 


P = I- 



I = 1% jkl = r (6*8 ji + fiilfijk) 


(89a) 
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(89b) 


Q = Qijki = -(DikSji + DaSjk + Dj k 8u + Djidik) — 2 DijD k i 


R = ZD ® D — (D ® 6 4- 8 0 D) 4- —6 <8)6 ( 89c) 

o 

with D = d®d\ d =unit vector defining the ’’preferred direction” or normal to the plane 
of isotropy and where a,/3,\,Gi, and Gt are elastic moduli and the symbol <g> represents 
a tensor product. 

Remark 3 For many materials, including metals and polymer-based compounds, per- 
haps the most physically-appealing separation of relaxation mechanisms, according to eqn. 
(88), is that for viscoelastic deviatoric and purely-elastic volumetric response; i.e., with 
finite and infinite relaxation times, respectively (r and p in separate ”integral” equations 
for each in eqns. 86 and 87). 

Remark 4 With a mew toward capturing a wider range of temperature a further gen- 
eralization, providing flexibility for delineating a time dependent and time independent 
regime (see Fig. 1 ), would be to explicitly take 

n* = ( - sun) (so) 

where the cut off value is taken as 

too(T) = 

The resulting implication would then be reflected by a distinct change in material modulus 
within the now partitioned quasilinear regime as illustrated in Fig. 6 and eqn. (85). 

Remark 5 It is interesting to note that the same format of additive-decomposition 
of ”equilibrium-plus-nonequilibrium” stresses and associated energy contributions (Eqns. 
19) will remain, even under conditions of large viscoelastic deformations (using second 
Piola Kirchhoff stresses). This is in conformity with the mathematically consistent and 
(physically-motivated) approach of multiplicative split of deformations, that is the de- 
formation gradient tensors (and related ” right” Green deformation tensors). Note that 
under conditions of small deformations, the latter (multiplicative) decomposition leads 
directly to its present "linearized ” form of eqn. 8. For example, see [23] for applica- 
tions to viscoelasticity; and [ 24 ] for alternatives of strain partitioning under more general 
conditions of large deformation elastoplasticity. 
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5.3 Thermodynamic Admissibility; Work and Dissipation In- 
equalities 

When general forms of hereditary convolution representations are utilized, the question 
of thermodynamic admissibility conditions becomes crucial. This subject has a long 
standing history in the literature on viscoelasticity [e.g., [8], [9], [25]-[28]]. For example, 
considering the one-dimensional case for simplicity (trivially generalized to the multi- 
axial case), Rabotnov [9] and Breuer and Onat [25], have formulated the sufficient (but 
not necessary) condition of positive work; wherein 


cr(t) = f ^(t- s) A f(s)ds 

Jo 

W = j a(s)‘j(s)ds > 0 

Jo 


(91) 

(92) 


and where a and 7 are the one-component (non-equilibrium) stress (see eqn. (87)) and 
strain field, respectively, and \&(t) is a general relaxation kernel. Further conditions 
involving the nonnegative definiteness and strict monotonicity and convexity of 'l'(t) 
have been also added: 


¥(s) > 0, tf'(s)<0, (93) 

\I>"(s) > 0, for any s > 0 (94) 

These are typically derived by considering step-loading strain functions. However, in 
order to define a complete set of conditions on $ (t) a more comprehensive set of conceived 
loading programs should be utilized; i.e., a sequence of many constant-strain steps of 
alternating signs, occurring at times t = 0, h, 2h, ...,nh in the limit h — * 0, thus leading 
to the use of the theory of distributions [[29], [30]] for generalized/singularity loading 
functions (delta, dipole, ... functions) 7 n (s) = <5 (n) = g (see Fig. 7 for an illustration). 
Together with the requirement of fading memory, it can be shown that this leads to 
more conditions of the type 

(_l) n tfl y( — > Q> s > 0, n = 0,l,2... (95) 

ds n 

The above condition is simply derived as follows. With 7 = $ n \the first of eqn. (92) 
gives a(t) = ^ n \t), whose sign must then be matched to that dictated by the fading- 
memory assumption. To this end, and depending on whether ”n is even (or odd), we 
note that a strain-rate history 7(3) = <5 (n) (s) is expected to give, for t > nh (in the limit 
h — »• 0), a no nn egative (nonpositive) stress response, in accordance with a fading memory 
hypothesis for a strain-rate history having identical (in magnitude, but opposite signs) 
changes, with the more recent being positive (or negative). In their present form, the 
above conditions (eqns. (95)) define a completely monotonic function. Although this 
excludes some of the nonstandard relaxation functions (e.g., used in instability studies of 
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shockwave fronts in viscoelastic composites [[28], [31]]), the above conditions will eliminate 
any physically unrealistic representations, such as "negative” viscosity in the equivalent 
mechanical analog of series of springs/dashpots [27]. It is important to note that utilizing 
the internal-variable formalism, i.e., with the complementary energy and dissipation rate 
functions discussed herein, all of the above conditions are automatically satisfied once 
the positive-definiteness and non-negativeness of the elastic ( Eij ki , My W ) and viscosity 
{Vijki) moduli are ensured, respectively. 

As an alternative to the above derivation, a more rigorous mathematical approach 
can be utilized here. Its formal setting is that for functional analysis of the problem of 
evolution with linear operators using the semigroup theory (e.g., [32]). For the present 
case, the semigroup corresponds to the convolution operator, T t+s = T t oT, for mapping 
composition, with {7V y — (^*7)(t)}, t > 0, where indicates convolution, and subscript 
find argument ”t” indicate the parameter for this one-parameter family {} of continuous 
linear operators. Aside from a number of technical conditions that are standard in the 
literature of semigroup theory, there are two important conditions here. The first stems 
from the existence of the natural, nonnegative, semi-norm induced by the associated 
inner product; i.e., here this corresponds to the positive work condition in eqn. (92) 

W =< T t 7,7 >=< <r ,7 >= f a(s)‘y(s)ds > 0 (96) 

Jo 

Using the history 7(s) = «5^ n ^(s), find carefully accounting for the singularity of 6^(s) 
in calculating the inner product, we use ( a(s) = \f^ n )) 

W =< a,7 >= (tf (n) *7)(0) 


or 

W = (*< n > * # n >)(0) = (-l) 2n ^ (2n) (0) > 0 (97) 

where 6(s) — 6(— s) = <5(s) = reflexive mapping of Dirac’s distribution or 6— function, 
and (0) indicates functions evaluated at t=0. Generalization from initial (zero) to any 
fixed (time a) time is trivial, through the convolution translation rule (^ * <5)(s — a) = 
('J'(s) * S(s — a)). Thus we get from (97), for any general time t, 

(_!)2"*(2n)(f) > 0 t > o (98) 

Following a similar procedure, the second condition, pertaining to the dissipativity 
(contractivity) of the associated in fini tesimal generator [32], here reduced simply to the 
dissipation inequality, e.g., see fly, in eqn. (9): 

Dy = - < 7, D t T t i > 

= -(7*<t)(0) 

= -(£ (n) * * 6 (n) )( 0) 
at 
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(99) 


or using convolution-time-shift properties 

= -$r( 2 " +1 >(*) > 0 

for any t > 0. 

The above conditions in (98) and (99) were obtained under very general conditions; 
they are applicable irrespective of the extent or type of material memory implied by the 
use of general kernels in eqns. (86) and (87). Remarkably, they have combined to give 
exactly the same conditions of the nonnegative (or nonpositive) characteristics for even 
(or odd) order time derivatives of the kernel functions. 

Remark 6 Starting from the additive decomposition for stress in eqn. (26), one 
can then adopt any general kernel in the equation defining the integral form of the 
nonequilibrium stress, cr m ( see eqn. ( 87) ), and proceed to use the constraint eqns. 
(98) and (99), with cr m and e replacing a and 8, respectively, as formal definitions 
for the stored energy (the complementary Gibb ’s function follows through a con- 
ventional Legendre transformation) and viscous dissipation function, respectively. 
This has two main advantages. Firstly, the direct definition, specially for the stored 
energy, will alleviate the well-known difficulties that may arise in attempting to de- 
rive them using alternative approaches (often with further assumptions), e.g., see 
[[8], [9]]. Secondly, satisfying the shown inequalities in eqns. (98) and (99) will also 
automatically ensure thermodynamic admissibility. 

Remark 7 In conjunction with the previous remark, and generalizing to multiaxial 
cases, i.e., with e(s) - <5 (n) (s)e 0 where e a is the fixed ” strain-direction ” tensors. We 
further need the positive (semi) definiteness of tensors E, M, rj (in addition to eqns. 
(98) and (99)) for full compliance with the thermodynamic admissibility conditions. 

6 Recursive Algorithm For Numerical Integration 

Considering the internal variables, the key aspect in the formulation of discrete time 
stepping procedures concerns the evaluation of the convolution integral in eqn. (86). In 
the present context, and utilizing the strain-driven form in eqn.(87), as typical in finite- 
element computations, the crucial observation is that the following recursive relation 
holds, for advancing the solution from step t n to t n +i with A t = tn+i — In '■ 
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tr™ +1 = e » cr™ + e s * rf Ae (100) 

where Ae = e„ + i — e n is the given strain increment in the step. In the second term of 
eqn.(100), the midpoint rule, known for its second-order accuracy! [33], [34]] has been 
used to linearly approximate the derivative. 

Note that the above algorithm requires the values of O’JJ* as the only historical data 
to be stored. Also, that this scheme provides for unconditionally stable[[33],[34],[35]] 
time stepping as well as the correct limits for large and vanishingly small times. 


7 Parametric Study 


Here we will utilize eqns. (54) and (56) to conduct creep parametric studies to study 
the influence of the various material parameters with an eye toward providing guidance 
and insight into the experimented characterization of the proposed model. Basically six 
material parameters will be assessed, i.e., E a ,El m ,E i/ m , and 1 /^; however, in keeping 
with our previous discussion (see eqns. (86) and (87)) E m and will be studied in the 
context of their respective material time clocks (i.e., relaxation and creep (retardation) 
times) pi = M -1 »7 and rl = EM 1 ff, respectively, and their ratio qI = (jj — 1)1 = E -1 M. 
Such a reduction can be accomplished only if we asstime all Poisson’s ratios to be equal, 
i.e., v = v s = v m = i/jj, and then take the ratio of the various moduli tensors as indicated 
above. In this way, eqns. ( 69) and ( 70) are reduced, such that in the case of uniaxial 
creep X and Y become: 


X = ff*(l - (^)e-J) 


Y = ^(l- (^)e-i) 

E, K T ’ ’ 


1 

0 

0 

1 

—v 

—V 


( 101 ) 


( 102 ) 


and the parametric study is reduced to the investigation of only four independent para- 
meters, that is, E a (or E),r,p, and v. In all creep test simulations the load (stress) was 
applied along the 11 direction, at a rate of 0.5 ksi/sec, to a constant level of 1 ksi and than 
held fixed for a given duration of typically lOr. The parametric study begins with the 
influence of one’s choice for the various Poisson’s ratios and is then followed by the effect 
of the infinite elastic stif fn ess (E 4 ) , the influence of the stif fn ess of the dashpot (viscous) 
element relative to the nonviscous stiffnesses (i.e., r), and then lastly the influence of the 
viscous to the nonviscous stiffness within the Maxwell element itself (i.e., p). Note that 
the variations of these parameters can be interpreted as simulating the behavior of the 
model over a given temperature range, as this would be a mechanism that could cause 
such a perturbation in the various moduli and Poisson’s ratios. 
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7.1 Effect of Poisson’s Ratios 

It is important to note that even with the restrictions to isotropy and uniaxial loading 
there remain in the general forms of eqns. (70) significant triaxiality and interaction 
effects among the various internal stress and strain components for nonzero and distinct 
Poisson’s ratios of the different elastic and viscous elements. To investigate this interac- 
tion behavior, we examined the seven cases given in Table 1, where the stiffness values 
were fixed at base line values (assumed based on Fig. 2b of Part II [10]) while the values 
of the Poisson’s ratios were changed. 


Table 1 Material parameters med for results generated in Figa. 8-18. 


Case 

E s (ksi) 

E m (ksi) 

p (sec) 

E„ (ksi-sec) 

r (sec) 

v a 

I'm 

v V 

1 

8500. 

8500. 

333.35 

2833475. 

666.7 

0 

0 

0 

2 

33 

33 

33 

33 

33 

0.2 

0.2 

0.2 

3 

33 

33 

33 

33 

33 

0.3 

0.3 

0.3 

4 

33 

3 ) 

33 

33 

33 

0.4 

0.4 

0.4 

5 

33 

33 

33 

33 

33 

0.49 

0.3 

0.3 

6 

33 

33 

33 

33 

33 

0.3 

0.3 

0.49 

7 

33 

33 

33 

33 

33 

0.3 

0.49 

0.49 


The creep results for cases 1 through 4 are shown in Figs. 8-12, where Fig. 8 shows the 
stress versus strain, Figs. 9 and 10 the 11- and 22-strain versus time history, Figs. 11 
the eq uili brium stress (stress in the spring element) in the 11 direction versus time, and 
finally Fig. 12 the effective overall Poisson’s ratio versus time. Examining these figures 
and eqns. (101) and (102), it is clear that when the individual Poisson’s ratios are all 
equal: 

• The stress versus strain behavior is unaffected, irrespective of the value of Poisson’s 
ratio taken, as is the strain in the 11 direction and the equilibrium stress in all 
directions. 

• The 22 (shown here) and 33 strain component is, as expected, highly influenced by 
the value of the Poisson’s ratio chosen. 

• Only primary creep response is generated, with this transient period ending after 
approximately 7 r, when the equilibrium stress in the loading direction becomes 
equal to the applied stress in that direction, i.e., the elastic limit of the material is 
reached. 

• The eq uili bri um stress in the transverse directions are zero, as expected given that 
the material is homogenous and the applied stress in those directions are zero. 

• The effective Poisson’s ratio is independent of time and equal to the Poisson’s 
ratios of the individual elements. 
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Considering three additional cases (i.e., cases 5, 6, and 7, respectively) in which, 1) 
the single spring element is taken to be purely deviatoric, 2) the viscous element in the 
Maxwell element is deviatoric and 3) the entire Maxwell element is taken to be purely 
deviatoric (i.e., v m — « 0.5). Note that in these cases the material is, in essence, 

implicitly assumed to be inhomogeneous. Consequently, the specialized uniaxial forms 
in eqns. (101) and (102) are not applicable and the more complex but general eqns.(69) 
and (70) must be used. The results for cases 1, 5, 6 and 7 of Table 1 are shown in Figs. 
13 through 18. Examining these figures one immediately sees that if the Poisson’s ratios 
are not fill taken equal to one another: 

• The initial effective stiffness can be increased (see Fig. 13). 

• The time required to reach the elastic limit can increase to as much as lOOr, as 
there is likely to be more than a single internal clock present (see G* and K*). 

• Although the total accumulated 11 strain, and 11 and 22 equilibrium stress com- 
ponents are the same as before, the transient portions are highly influenced by the 
assumed value of the various Poisson’s ratios. In particular note the difference in 
the sign of the equilibrium stress in the 22 direction for case 5 as compared with 
cases 6 and 7. (see Figs. 14, 16 and 17) 

• The effective Poisson’s ratio becomes time dependent for the cases when the indi- 
vidual element Poisson’s ratios are not the same. With the initial effective ratio 
being equal to that given in eqn. (78) and the final being that of the single spring 
element (i.e., i/ a ). This observation provides us with insight into the needed critical 
test to characterize, as well as validate, the present model construction, (see Fig. 
18) 

Obviously, this demonstrates the rather limited usefulness of one-dimensional mod- 
eling and fitting in realistic representations for general applications. In fact, any ad-hoc 
extrapolation assumptions made [e.g., [11] or [36]] can lead to rather erratic history pre- 
dictions of the implied Poisson’s-effect history (i.e., time intervals of alternating negative 
and positive Poisson’s ratios). Apparently the difficulty stems from the complex depen- 
dency of both the transverse as well as the axial strain components on the three coupled 
Poisson’s factors, thus rendering it extremely difficult to define a simple operator form 
(creep compliance/ relaxation modulus) for their ratios as would be the case in these 
simple ad-hoc extrapolation schemes. Such difficulties and anomalous predictions can of 
course be completely by-passed provided the coaxiality of all elastic moduli tensors are 
imposed and the conventional engineering material parameter bounds, i.e., positive mod- 
uli and 0 < Vi < 0.5 are adhered to in our general formulation (ensuring the convexity 
of the two complementary potentials under study). This coaxiality assumption appears 
to be experimentally justifiable at this time, as preliminary measurements (see Part II, 
[10]) have indicated that the effective Poisson’s ratio is time independent. 
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7.2 Effect of Infinite Elastic Stiffness, E s 

Next, four cases were examined in which the infinite elastic stiffness (i.e., that of the lone 
spring element, in Fig. 4a) was varied from 2000 to 12500 ksi as identified in Table 2. 
Note that in these cases, as in the previous seven, the spring stiffness (E a ) is taken equal 
to that of the spring stiffness within the Maxwell element (i.e., Em), that is the ratio 
a — 1. Results are presented in Figs. 19, and 20, where the resulting stress versus strain 
history and strain versus 


Table 2 Material parameters used for results in Figs. 19-20. 


Case 

E a (ksi) 

Em (ksi) 

p (sec) 


r (sec) 




1 

2000. 

2000. 

333.35 

666700. 

666.7 

0.3 

0.3 

0.3 

2 

4500 

4500 

)> 

1500075. 


0.3 

0.3 

0.3 

3 

8500 

8500 


2833475. 


0.3 

0.3 

0.3 

4 

12500 

12500 


4166875. 


0.3 

0.3 

0.3 


time history are shown. Figures 19 and 20 clearly indicate, as one might expect from 
eqn. (102), that the effective load-up stiffness is strongly influenced (increased) and the 
accumulation of time-dependent strain is also highly influenced (decreased) by increasing 
the assumed value of E s . Such a result is necessitated by our insistence upon maintaining 
a constant p, r and a ratio. Similarly, the resulting equilibrium stress versus time history 
and effective Poisson’s ratio history are identical with those described in Figs. 11 and 12, 
respectively, and the total accumulated transverse strains are clearly just the opposite of 
those shown in Fig. 20 with the magnitude being scaled by the Poisson’s ratio (see eqn. 
(102)). As an aside, the trends observed in Figs. 19 and 20 are in keeping with those 
one might expect as they decreased the temperature of the material. 


7.3 Effect of Varying the Characteristic Creep (Retardation) 
Time, r 

Now we consider four cases which illustrate the influence of r on the creep response of 
the material, wherein all other parameters are fixed at the base line values. As r is a 
ratio (EM _1 rj) and we desire that p remain constant, E*,, must therefore be modified. 
The four sets of material parameters used are given in Table 3. Note that by necessity a 
is equal to 0.1, 1, 2, and 5 for cases 1, 2, 3, and 4, respectively. 


TVK)e 3 Material par^wigten naed for results in Figs. 21-22. 


Case 

E 3 (ksi) 

Em (ksi) 

p (sec) 

E„ (ksi-sec) 

r (sec) 

v» 

V m 


1 

8500. 

850. 

333.35 

283347.5 

366.7 

0.3 

0.3 

0.3 

2 

8500 

8500 

r> 

2833475. 

666.7 

0.3 

0.3 

0.3 

3 

8500 

17000 

n 

5666350. 

1000. 

0.3 

0.3 

0.3 

4 

8500 

42500 

w 

14167375. 

2000. 

0.3 

0.3 

0.3 
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The resulting stress-strain and strain-time (and equilibrium stress (X, a a ) versus time) 
response histories are shown in Figs. 21 and 22, respectively. Examining eqn. (102) 
and Fig. 21, it is clear that the effective load-up modulus is once again increased as Em 
(or correspondingly a) and r are increased, however, this time the total accumulated 
strain remains the same in all cases; as the elastic limit strain (gj) remains fixed, i.e., 
E 5 = 8500 ksi. This is opposite to the cases in Fig. 19, where E, is continuously being 
increased and thus the elastic limit strain is being reduced in each case. Consequently, 
one has sufficient freedom in this model to modify the effective load-up modulus without 
correspondingly reducing the total amount of strain accumulation. Figure 22 and eqn. 
(102) again, clearly indicate another influence of increasing r, that is, to increase the 
duration of the primary (transient) zone of the time-dependent straining. Note that in all 
cases the primary zone is complete (with respect to the external clock) by approximately 
7 times the associated internal clock (i.e., r). This fact can be extremely helpful in 
characterizing the proposed model, in that we can deduce the internal clock from the 
experimentally observed saturation time (external clock). Similarly, if we experimentally 
observe that the saturation time of primary creep is independent of temperature, one 
would conclude that r should also be taken to be temperature independent. Lastly, as the 
individual Poisson’s ratios are all equal, the effective Poisson’s ratio is time independent 
(see Fig. 12) and the 22 and 33 components of total strain are simply appropriately 
modified mirror images of the 11 component, and the 22 and 33 equilibrium stresses are 
zero (see eqn. (102)). 

7.4 Effect of Varying the Characteristic Relaxation Time, p 

Finally, we will consider four cases which illustrate the influence of p on the creep response 
of the material, wherein all other parameters are fixed at the base line values. Once again, 
as p is a ratio (M _1 »7) and we desire that r and a remain the same as in the previous 
cases, i.e., r = 666.7 and a = 0.1, 1, 2, and 5 ; 17, must therefore be modified. The four 
sets of material parameters used to generate Fig. 23 are given in Table 4. 


Table 4 Material parameters naed for results in Figs. 23-26. 


Case 

E, (ksi) 

E m (ksi) 

p (sec) 

E,, (ksi-sec) 

t (sec) 

^5 

Vm 

El 

i 

8500. 


606.09 

515176.5 

666.7 

0.3 

0.3 

0.3 

2 

8500 

8500 

333.35 

2833475. 

n 

0.3 

0.3 

0.3 

3 

8500 

17000 

222.23 

3777910. 

5 > 

0.3 

0.3 

0.3 

4 

8500 

42500 

111.12 

4722600. 


0.3 

0.3 

0.3 


Here the resulting stress-strain history is identical to that of Fig. 21, as both E, and 
Em are the same as in the previous cases. The strain-time and equilibrium stress versus 
time histories are influenced however, as one might suspect, for example see Fig. 23. 
Comparing Figs. 23 to 22, we see no influence of p on the duration of the primary creep 
regime, as one would expect given eqn. (102). However, the rate of accumulation is 
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slightly influenced, given the same value of a, as seen in Fig. 23. Clearly, in relaxation 
p would now have a significant influence on the duration of stress relaxation. 

Revisiting the influence of the individual Poisson’s ratios, we re-examine the cases in 
Table 4, but this time we take v s = 0.1 while v m == u v = 0.3. Again, the stress-strain 
history is identical and the 11 component histories of the total strain and equilibrium 
stress are very si mil ar, the duration of the transient zone just being extended. The 22 
components, however, become highly time dependent as does the effective Poisson’s ratio 
of the material model. These response histories are illustrated in Figs. 24-26, respectively. 
Note that the behavior of the components transverse to the load are quite complex, and 
again the only sure way to validate the model assumptions is to collect experimental data 
in those directions. 


8 Conclusions 

Traditionally, the total strain has been partitioned into elastic (time-independent re- 
versible), inelastic (irreversible) and thermal (reversible) strain components, where the in- 
elastic strain corresponds to such physical phenomena as time-independent plastic strain 
or time-dependent viscoplastic strain (sometimes referred to as creep strain). In this 
paper alternative forms for the hereditary integrated representations of the viscoelastic 
response were discussed to account for the time/rate dependent reversibility typically 
observed in experiments on materials exhibiting quasilinear domains. This discussion 
was in the context of a general framework for the viscoelastoplastic material modeling 
of the internal-state- variable type. It was shown that such a potential based framework 
necessitated the specific partitioning of both the strain (reversible and irreversible) and 
stress (equilibrium and non-equilibrium) state variables. Important key issues, pertaining 
to both the theoretical foundation (e.g., thermodynamic admissibility) as well as the nu- 
merical implementation (e.g., recursive integral forms) of these general descriptions were 
discussed. In particular, explicit constraints on the underlying kernel functions were 
given. Furthermore, even though numerous multiaxial formulations have been noted in 
the literature in their symbolic form, it was demonstrated here that without the assump- 
tion of equal Poisson’s ratios (coaxiality of the various moduli tensors) the uniaxial 
reduction to the classical linear solid element, is not possible given the starting view- 
point of internal state variables, defined through either the differential or convolution 
integral approach. Of course, as implied by classical forms, taking all Poisson’s ratios to 
be zero is merely a special case of this equality condition. 

Several parametric studies and qualitative response assessments were made to help 
identify important factors in the actual characterization of the specific-form of the model 
studied for the titanium alloy, TIMETAL 21S, as discussed in Part II [10]. The more 
important, with respect to characterization, parametric results are itemized as follows: 

• If the experimentally observed effective Poisson’s ratio is independent of time 
then all individual mechanical elements’ Poisson’s ratios can be assumed equal 
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to that of the measured one, thereby significantly reducing the complexity of the 
required reversible constitutive model. Several other considerations (e.g., damping- 
extent ordering) were also found to favorably support this assumption. 

• Only primary creep response is generated, with this transient period ending after 
approximately 7 r (given the assumption of equal Poisson’s ratios), when the 
equilibrium stress in the loading direction becomes equal to the applied stress in 
that direction, i.e., the elastic limit of the material is reached. Consequently, r , the 
internal material clock, can be immediately obtained by dividing the experimentally 
observed saturation time by seven. However, this is not the case when different 
Poisson’s ratios are assumed wherein the period now intricately depends upon these 
ratios (and the ’’implied” number of internal clocks are increased). 

• The effective load-up stiffness is strongly influenced (increased) and the accumu- 
lation of time-dependent strain is also highly influenced (decreased) by increasing 
the assumed value of E a . Consequently, variation with temperature of E s should 
be a key factor under nonisothermal conditions. 

• If the assumption of constant Poisson’s ratios is accurate, than eqn. (79) proves 
that the dynamic modulus is equal to the sum of E„ and E m , i.e., Eq = E 3 + E m , 
and E m can be deduced from directly measurable quantities. 

• One has sufficient freedom in this model to modify the effective load-up modulus 
without correspondingly reducing the total amount of strain accumulation. Increas- 
ing r increases the duration of the primary (transient) zone of the time-dependent 
straining. Similarly, if we experimentally observe that the saturation time of pri- 
mary creep is independent of temperature, this would than imply that r should 
also be taken to be temperature independent. 

Armed with these insights, and assuming the validating of a time-independent ef- 
fective Poisson’s ratio the characterization of the. proposed multiaxial reversible model 
is straight forward. Alternatively, if the effective Poisson’s ratios are time dependent, 
then one may need the assistance of an automated material parameter estimator, [37], 
as would be necessarily the case for a wider relaxation mechanism spectrum (greater 
frequency-dependency of the moduli) calling for several internal state variables. The 
companion paper (see Part II [10]) will focus on the actual characterization and compar- 
ison with test results for a model material, i.e., TIMETAL 21S, over a wide range of 
temperatures (i.e., 23 to 650 °C). 

Finally, given the host of viscous effects present in structural applications at elevated 
temperatures, the practical utility of such a generalized hereditary form as presented 
herein is significant. For example considering that most aerospace engine designs are 
typically limited to the quasilinear stress and strain regimes, as in this case the reversible 
time-dependent response component becomes dominate in comparison to the irreversible 
component. Alternatively, one can envision another extreme case (e.g., in polymer and 
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rubber based systems under varying temperatures) in which a purely reversible viscous 
response is present, thus necessitating the use of such a theory. And lastly, an obvious 
natural extension for general applicability is the middle ground in which a combined 
reversible and irreversible representation is required. Here the inclusion of a reversible 
time-dependent response as that discussed in this paper will extend the applicable pre- 
dictive stress and/or temperature window of application. 
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(c) Low temperature. 



Time independent 
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Time dependent 
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Fig 2. — Rate sensitivity of static stiffness, timetal 21 S. 



(a) Generalized four element: viscoplastic. 


Partitioning the strain 
e = e e + c c 



€ C = € P + € S 


• P _ (g-«) n 
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(b) Classic four element: viscoelastic. 


Fig 1 . — Reversible and irreversible Fig 3.— Four element mechanical models depicting 

threshold surfaces. the deformation response of a material. 
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Partitioning the stress 



(a) Generalized standard solid. 
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(b) Solution response. 

Fig 4. — Three element mechanical model and corresponding creep and 
relaxation response history. 


e = e R + e ,R 



Fig 5. — General hereditary behavior model indicated 
with appropriate partitioning of stress and strain. 
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Irreversible 



e R = E' 1 <t t ,(T) + E” 1 <cr-cr T |(T)^> 

Fig 6. — Uniaxial depiction of reversible and irreversible 
domains. 
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Fig 7. — Step-strain loading approximating functions 

7 (n) = 2 H) n - k Y(kh) in the limit h-o. 

h n k =c k!(n-k)! 

According to fading-memory assumption, stress 
relaxes to positive or negative values for even (Case 
a, n = 2) or odd (Case b, n = 3) n- th order derivatives. 


Strain, y(t) 




(b) Case n = 3 (odd). 
Fig 7. — Concluded. 
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Fig 8. — Stress versus strain histories in the 1 direction 
given equal Poisson's ratios, see cases 1 -4 of Table 1 . 



Time, sec 


Fig 9. — Total reversible strain in the 1 direction versus 
time given equal Poisson's ratios, see cases 1 -4 of 
Table 1 . 



Time, sec 


Fig 10. — Total reversible strain in the 2 and 3 direction 
versus time given equal Poisson's ratios, see cases 
1 -4 of Table 1 . 



Fig 11.— Equilibrium stress (X) in the 1 direction versus 
time given equal Poisson's ratios, see cases 1 -4 of 
Table 1. 



Time, sec 


Fig 12. — Effective overall Poisson's ratio versus time 
given equal Poisson's ratios, see cases 1 -4 of 
Table 1. 



Fig 13. — Applied stress versus strain histories in the 
1 direction given nonequal Poisson's ratios, see 
cases 1 , 5, 6 and 7 of Table 1 . 
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T otal strain, Y 2 , Y 3 , (xl 0" 6 ) Total strain, Y, , (xl 0“®) Total strain, Y-, , (xl 0“®) 



0 1000 2000 3000 4000 5000 6000 7000 

Time, sec 



0 20,000 40,000 60,000 


Time, sec 

Fig 14. — Total reversible strain in the 1 direction 
versus time given nonequal Poisson's ratios, see 
cases 1 , 5, 6 and 7 of Table 1 . (a) 10 7. (b) 100 7. 



0 1000 2000 3000 4000 5000 6000 7000 


Time, sec 

Fig 15. — Total strain in the 2 and 3 directions versus 
time given nonequal Poisson's ratios, see cases 1 , 
5, 6 and 7 of Table 1. (a) 10 t. (b) 100 t. 



Time, sec 


Fig 15. — Concluded, (b) 100 t. 



Time, sec 

Fig 16. — Equilibrium stress in the 1 direction versus 
time given nonequal Poisson's ratios, see cases 1 , 
5, 6 and 7 of Table 1. (a) 10 t. (b) 100 t. 



Time, sec 

Fig 16. — Concluded, (b) 100 t. 
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Spring stress, X 2 , X 3 , ksi Spring stress, X 2 , 



Time, sec 


ig 17. — Equilibrium stress in the 2 and 3 directions 
versus time given nonequal Poisson's ratios, see 
cases 1 , 5, 6 and 7 of Table 1 . (a) 1 0 t limit. 

(b) 100 t limit. 



Time, sec 


Fig 18. — Effective overall Poisson's ratio (V 12 ) versus 
time given nonequal Poisson's ratios, see cases 1 , 
5, 6 and 7 of Table 1 . (a) 1 0 t limit, (b) 1 00 t limit. 



Fig 17. — Concluded, (b) 100 t limit. 
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Fig 19. — Influence of E s on stress (L) versus strain (Y) 
in the 1 direction, see Table 2. 



Time, sec 


Fig 20. — Influence of E s on total strain (Y) in the 1 
direction verses time, see Table 2. 



Fig 21 . — Influence of t on stress (L) versus strain (Y) 
history in the 1 direction, see Table 2. 



Time, sec 


Fig 22. — Influence of t on total strain (Y) in time 
response, see Table 3. 



Fig 23. — Influence of p on total strain (Y) in time 
response, see Table 4. 



Fig 24. — Influence of p on total strain histories in the 
2 and 3 directions given different Poisson's ratios 
in the spring and Maxwell elements. 
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Time, sec 

Fig 25.— influence of p on equilibrium stress (X) 
versus time in the 2 and 3 directions given different 
Poisson's ratios in the spring and Maxwell elements. 



Fig 26.— Influence of p, given different Poisson's ratios 
in the spring and Maxwell elements, on the effective 
Poisson's ratio (V 12 ) history. 
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